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1. Introduction 

The accurate and efficient representation of shapes is a major challenge in geometric 
modeling. To achieve high order accuracy in the representation of curves, surfaces 
or functions, piecewise polynomials models are usually employed. Parametric models 
with prescribed regularity properties are nowadays commonly used in Computer Aided 
Geometric Design (CAGD) to address these problems. They involve so-called spline 
functions, which are piecewise polynomial functions on intervals of M with continuity 
and differentiability constraints at some nodes. Extensions of these functions to higher 
dimension are usually done by taking tensor product spline basis functions. Curves, 
surfaces or volumes are represented as the image of parametric functions expressed in 
terms of spline basis functions. For instance, surface patches are described as the image 
of a piecewise polynomial (or rational) map from a rectangular domain of to M”. 
But to represent objects with complex topology, such maps on rectangular parameter 
domains are not sufficient. One solution which is commonly used in Computer-Aided 
Design (CAD) is to trim the B-spline rectangular patches and to “stitch” together the 
trimmed pieces to create the complete shape representation. This results in complex 
models, which are not simple to use and to modify, since structural rigidity conditions 
cannot easily be imposed along the trimming curve between two trimmed patches. 

To allow flexibility in the representation of shapes with complex topology, another 
technique called geometric continuity has been studied. Rectangular parametric surface 
patches are glued along their common boundary, with continuity constraints on the 
tangent planes (or on higher osculating spaces). In this way, smooth surfaces can 
be generated from quadrilateral meshes by gluing several simple parametric surfaces, 
forming surfaces with the expected smoothness property along the edges. 

This approach builds on the theory on differential manifolds, in works such as [7], 
[12], [n|- The idea of using transition maps or reparameterizations in connection with 
building smooth surfaces had been used for instance by DeRose [7| in CAGD, who gave 
one of the hrst general dehnitions of splines based on hxing a parametrization. 

Since these initial developments, several works focused on the construction of such 
surfaces [IH], [IE], |2l], [23], [5], [32], [H], [13], [lO], [27], [26], [3], ... with polyno¬ 
mial, piecewise polynomial, rational or special functions and on their use in geometric 
modeling applications such as surface htting or surface reconstruction 0 , 1251 , 119 ,... 

The problem of investigating the minimal degree of polynomial pieces has also been 
considered [19]. Other research investigates the construction of adapted rational tran¬ 
sition maps for a given topological structure [T]. We refer to [20] for a review of these 
constructions. Constraints that the transition maps must satisfy in order to dehne 
regular spline spaces have also been identihed [21] • But it has not yet been proved that 
these constraints are sufficient for the constructions. 

The use of spline functions to approximate functions over computational do¬ 
mains with arbitrary topology received recently a new attention for applications in 
isogeometric analysis. In this context, describing the space of functions, its dimension 
and adapted bases is of particular importance. A family of bi-cubic spline functions 
was recently introduced by Wu et al [31] for isogeometric applications, where constant 
transition maps are used, which induce singular spline basis functions at extraordinary 
vertices. Multi-patch representations of computational domains are also used in [3], 
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with constant transition maps at the shared edges of rectangnlar faces, using an iden- 
tihcation of Locally Rehned spline basis functions. In [13], continuous splines are 
described and the condition is transformed into a linear system of relations between 
the control coefficients. The case of two rectangular patches, which share an edge is 
analyzed experimentally. In [2] , the space of G^ splines of bi-degree ^ 4 for rectangu¬ 
lar decompositions of planar domains is analyzed. Minimal Determining Sets of points 
are studied, providing dimension formulae and dual basis for G^ spline functions over 
planar rectangular meshes with linear gluing transition maps. 

Our objective is to analyze the space of G^ spline functions for rectangular and 
triangular patches with arbitrary topology and general rational transition maps. We 
are interested in determining the dimension of the space of G^ spline functions which 
are of degree ^ k on triangular pieces and of bi-degree ^ {k, k) on rectangular pieces. 
To define the space of G^ spline functions, we introduce the concept of topological 
surface with gluing data attached to the edges shared by the faces. The framework 
does not require manifold constructions and is general enough to allow non-orientable 
surfaces. We describe compatibility conditions on the transition maps so that the 
space of differentiable functions is ample and show that these conditions are necessary 
and sufficient to construct ample spline spaces. A separability property is involved to 
obtain a dimension formula of the G^ spline spaces of degree ^ k on such topological 
surfaces, for k big enough. This leads to an explicit construction of basis functions 
attached respectively to vertices, edges and faces. 

For the presentation of these results, we structure the paper as follows. The next sec¬ 
tion introduces the notion of topological surface Ad, differentiable functions on Ad and 
constraints on the transition maps to have an ample space of differentiable functions. 
Section 3 deals with the space of spline functions which are piecewise polynomial and 
differentiable on Ad. Section 4 analyzes the gluing conditions along an edge. Section 
5 analyzes the gluing condition around a vertex. In Section 6 , we give the dimension 
formula for the space of spline functions of degree ^ k over a topological surface Ad 
and describe explicit basis constructions. Finally, in Section 7, we detail an example 
with boundary edges and another one with no boundary edges. We also provide an 
appendix with an algorithmic description of the basis construction. 

2. Differentiable functions on a topological surface 

Typically in CAGD, parametric patches are glued into surfaces by splines (i.e., poly¬ 
nomial maps) from polygons in The simplest G'^ construction is with the polygons 
in situated next to each other, so that G'^ continuity across patch edges comes 
from continuity of the coordinate functions across the polygon edges. This is called 
parametric continuity. A more general construction to generate a G^ surface from 
polygonal patches is called geometric continuity 0 , 123 . Inspired by differential ge¬ 
ometry, attempts have been made lEI, |2g, H to dehne geometrically continuous G^ 
surfaces from a collection of polygons in with additional data to glue their edges 
and differentiations. They are dehned by parametrization maps from the polygons to 
satisfying geometric regularity conditions along edges. 

It is easy to dehne a G^ surface from a collection of polygons and homeomorphisms 
between their edges. 
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Definition 2.1. Given a collection M 2 of (possibly coinciding) polygons cTj in a 
topological surface M is defined by giving a set of homeomorphisms fi : Ti ^ Tj between 
pairs of polygonal edges Ti C cTj, Tj C aj (cri,aj G M2)- Each polygonal edge can be 
paired with at most one other edge, and it cannot be glued with itself. 

A -continuous function on the topological surface M is defined by assigning a 
continuous function fi to each polygon (jj, such that the restrictions to the polygonal 
edges are compatible with the homeomorphisms pi. 

The topological surface M is the disjoint union of the polygons, with some points 
identified to equivalence classes by the homeomorphisms pi. The polygons are also 
called the faces of M and their set is denoted M 2 . Each homeomorphism p, identifies 
the edges r*, Tj of the polygons Uj, aj to an interior edge of M. We say that the edge 
is shared by the faces cxj and aj. An edge not involved in any homeomorphism /i is a 
boundary edge of M. The edges of M are the equivalent classes of edges of the polygons 
of identified by the homeomorphisms p. Their set is denoted Mi. Similarly, let 
Ml) denote the set of Ad-vertices, that is, equivalence classes of polygonal vertices. 
An interior vertex is an equivalence class of polygonal vertices 70 , 71 ,..., 7 ^ = 70 such 
that the adjacent vertices 7 j, 7 j+i are identified by an edge homeomorphism. The set 
of these equivalent classes of identified vertices of the polygons, or interior vertices of 
M, is denoted Mq. 


2.1. Gluing data. The definition of a differential surface S typically requires an atlas 
of S', that is a collection {VpMp}p&j such that {Vp}p^j is an open covering of S [50] . 
Each fjp is a homeomorphism ipp -.Up ^ Vp, where Up is an open set in For distinct 
p,q e J such that V), fl Eg 7^ 0, let Up^q := 'ipp^iVp fl Vq) and Uq^p := fl Vq). 

Then the map o 'ipp : Up^q —>■ Uq^p is required to be a G^-diffeomorphism. The 
maps cppq : 'ip~^ o ipp are called transition maps. A differentiable function / on S' is a 
function such that for any open set Vp, the composition fp = fo -0“^ : f/p C —)■ M 
is differentiable. 

Our objective is to study the space of differentiable functions that can be constructed 
on a surface S' associated to the topological surface M. Instead of an atlas of a 
differential surface S, we consider a topological surface M together with gluing data 
given by maps (that we call transition maps) between the pairs of faces of M that 
share an edge in M. We make this precise in the following definition. 


Definition 2 . 2 . For a topological surface M, a gluing structure associated to M 
consists of the following: 


A 


for each face a G M 2 an open set of containing cr; 

for each edge r G Adi of a cell a, an open set Ur,a of containing r; 

for each edge t E Mi shared by two faces ai,aj G Ad 2, a GEdiffeomorphism 

called the transition map '■ Ut,!^ — t Ur,aj between the open sets Ur^cn and 

Ur,(rp and also its correspondent inverse map 

for each edge r G Adi of a cell a, the identity G^-diffeomorphism that defines 

= Id between U„ and U^ 


the, identity transition .rpap,-^ 
transition map as m Dehmtion 2.2 


r,= id between and Ur.,a- . . 
ditters from the usual notion of transition map 

in the context of differential manifolds (see [ 12 ]), since we do not require compatibility 

conditions at the vertices and across edges. The precise compatibility conditions that 

we need on these maps (po 


are given in Sections 2.3 ans 2.4 
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Figure 1 . Topological surface constructed from two triangles. 


Let r = (Ti,r2) be an edge shared by two faces cti,ct2 G M2 and let 7 = (71,72) 
be a vertex of r corresponding to 71 in ai and to 72 in a2, as in Figure We denote 
by (resp. r^) the second edge of di (resp. CT2) through 71 (resp. 72) We associate 
to (Ti and (72 two coordinate systems (Mi,ni) and {u2,V2) such that 71 = (0,0), Ti = 
{(mi, 0 ),mi e [ 0 , 1 ]}, t[ = {( 0 ,ni),ni G [ 0 , 1 ]} and 72 = ( 0 , 0 ), T2 = {(0,^2), 1^2 G [ 0 , 1 ]}, 
T2 = {{u 2 , 0 ),U 2 G [ 0 , 1 ]}. Using the Taylor expansion at ( 0 , 0 ), a transition map from 
Ur,ai to Ur,a2 IS then of the form 


( 1 ) 


^ 2,^71 : {ui,Vi) —)■ (^^ 2 ,^^ 2 ) 


/ Vibr,-y{ui) + vfpi{ui,vi) \ 

[ui + Vi Ur,7(^1) + vIp 2 {Ui, Ui) J 


where ar,^{ui), br^-y{ui), pi{ui,Vi), p2{ui,Vi) are functions. We will refer to it as the 
canonical form of the transition map 0o-2,(ti 7 along r. The functions are 

called the gluing data at 7 along r on cxi. 


Definition 2.3. An edge t & M which contains the vertex 7 G Ai is called a crossing 
edge at 7 if aT-^.^( 0 ) = 0 where [a,-, b,-] is the gluing data at 7 along r. We dehne 
c.,-(7) = 1 if r is a crossing edge at 7 and c,-(7) = 0 otherwise. By convention, £,-(7) = 0 
for a boundary edge. If 7 G Aio is an interior vertex where all adjacent edges are 
crossing edges at 7, then it is called a crossing vertex. Similarly, we dehne C4.(7) = 1 
if 7 is a crossing vertex and c+(7) = 0 otherwise. 

2 . 2 . Differentiable functions on a topological surface. We can now dehne the 
notion of diherentiable function on M: 


Definition 2.4. A diherentiable function / on the topological surface Af is a collection 
/ = {fa)(7eM of diherentiable functions /„- : f/o- —)■ M such that V7 G r = ai fl (J2, 
VU G Ur,ar , 

(2) J^(/<,J(u) = o 0^2^^J(u) 

where J7 is the jet or Taylor expansion of order 1 at 7. 

If fi, /2 are the functions associated to the faces cxi, (T2 G Af 2 which are glued along 
the edge r with a transition map of the form ([^, the regularity condition ([^ leads to 
the following relations: 

• /i(mi, 0 ) = /2 o0a2,<7i(Mi,O) for Ui G [ 0 , 1 ]; that is 

/i(M1 ,0) = /2(0 ,Mi) 


( 3 ) 
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(4) 


dfi, . 

to 


dh 

dv\ 


d{f2 o 0) 

dvi 


(mi,0) = 


(mi,0) 




for 0 = 00-2,0-1 and Ui E [0,1], which translates 

df2.^ X , .X df2 , . 

dU2 dV2 


for Ml e [0,1], with a(Mi) = |^(mi, 0), b(Mi) = §^(mi, 0), where 0i and 02 
are the components of 0 at the hrst and the second variable respectively. 

A convenient way to describe this regnlarity condition is to express the relation Q as 
a relation between differentials acting on the space of differential fnnctions on the edge 
r: 


(5) ci,-,^('Ui)0i,2 A ^vi 0 

With this notation, at a crossing vertex 7 with 4 edges we have = 0- The 

differentials along two opposite edges are “aligned”, which explains the terminology of 
crossing vertex. 

Definition 2.5. A snbspace T> of the vector space of differentiable fnnctions on M. is 
said to be ample if at every point 7 of a face a of Ad, the space of valnes and differentials 
at 7, namely [/(y), 0„,(/)(7), d^Ji^)] for f eV, is of dimension 3. 

This dehnition does not depend on the choice of the face a to which 7 belongs, since 
for 7 on a shared edge, the valne and differentials coincide after transformation by the 
invertible transition map. 


2.3. Compatibility condition at a vertex. Giving glning data on the edges is not 
snfiicient to ensnre the existence of an ample space of differentiable fnnctions on Ad. 
At vertices shared by several edges and faces, additional conditions on the transition 
maps need to be satished. We describe them in this section, and show that they are 
snfiicient to constrnct an ample space of splines on Ad in the following sections. 

For a vertex 7 G Adg, (see Fig. which is common to faces (Ti,...,(Ti? glned 
cyclically aronnd 7, along the edges r* = (Tj+i fl (Tj for z = 1,..., F (with ap+i = Ui), 
we impose the following condition: 

(6) O ■ ■ ■ O J.^(03,2) O J.^(02,i)(m, v ) = (m, v ) , 

where Ay is the jet or Taylor expansion of order 1 at 7. We can assnme that for each 



Figure 2. The faces cij for z = 1,..., 5 are glned cyclically aronnd a 
vertex 7. 
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i = 1,..., F, the edge r, is defined (linearly) by = 0 in cxj. It is easy to check that 
the condition ([^ on the Taylor expansion at 7 leads to the following: 

Condition 2.6. If the vertex 7 is on the faces ai,... ,aF glued cyclically around 7 , 
the gluing data [a*, bj] at 7 on the edges Ti between (Ti_i and Oi satisfies 


( 7 ) 



1 

ai( 0 ) 



This gives algebraic restrictions on the values ai(0), bj(0). At a crossing vertex 7 
(see Def. 2.3) of four incident edges, the equality ([^ amounts to 


( 8 ) 


bi( 0 )b 3 ( 0 ) = l, 62 ( 0 ) 64 ( 0 ) = 1. 


It turns out that Condition |2.6| is not sufficient around crossing vertices for ensuring 
an ample space of differentiable functions on Ad. An obstruction was noticed in IZB 
in a setting of rectangular patches. We write this constrain in a general setting: 


Condition 2.7. If the vertex 'y is a crossing vertex with 4 edges ti, ... ,T 4 , the gluing 
data [oj, bj] z = 1... 4 on these edges at 7 satisfy 


( 9 ) 

( 10 ) 


a'l ( 0 ) + 

“ 2 ( 9 ) + 


q(o) 

64 ( 0 ) 

TM 

f.,(0) 


-b.(0) 

-b2(0) 


(“ 3 ( 9 ) + 

(^ 04 ( 9 ) + 


b' 2 ( 0 ) \ 

b 2 ( 0 )t 

b'3(0) \ 

bs(0)t 
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Lemma 2.8. If the space of differentiable functions on A4 is ample and Condition \2.6 
is satisfied, then the gluing data at every crossing vertex 'y of 4 incident edges must 
also satisfy Condition\2.7 


Proof. The value and first derivatives at every point 7 G Ad of all differentiable func¬ 
tions on Ad should span a space of dimension 3. 

If 7 is a crossing vertex, then we have 4 restrictions on the Taylor expansions of a 
spline components (/i, / 2 , /s, f^). Let us write the Taylor expansion of /, at 7 = (0, 0) 
as fi = Pi + QiUi + riVi -\- SiUiVi + .... The gluing conditions imply the following. From 

(§, 

(11) Pi=P2=P3= P4, qi = 7-2, q 2 = 7 - 3 , qs = and ^4 = t^i 

this together with the Condition (|^ on the first derivatives imply 

t’2(0)gi - g3 = 0, 63 ( 0 )g 2 - = 0. 

When we consider the derivative of fi with respect to Ui, again applying (|^, we get 
the conditions 


52 — 62(0)54 — 62(0)^! -|- a2(0)Q'2, 

53 — 63(0)82 = 63(0)g2 + ci3(0)q'3, 

54 - 64(0)53 = 64(0)g3 -h 04(0)^43 
Si - 64(0)54 = 6'4(0)g4 + ai( 0 ) 9 i- 
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By combining the last four equations respectively with the weights b 4 ( 0 ), bi( 0 )b 4 ( 0 ), 
bi(0), 1, together with Condition 2.6, we get 


( 12 ) 


Qi (cii(O) + b 2 ( 0 )b 4 ( 0 )) + ga bi( 0 ) (b 2 ( 0 )a 3 ( 0 ) + b 2 ( 0 )) + 

q 2 (bi(O)a'(0) + b;(0)) + q, b 2 ( 0 )(al( 0 ) + bi( 0 )b' 3 ( 0 )) = 0. 


This relation does not involve the cross derivatives si, S 2 , S 3 , S 4 , but gives an unwanted 
relation between the hrst order derivatives. After replacing ([IT| and ([^ in (12), we 


encounter conditions ([^ and (10). Under these conditions, there is no relation between 
qi,q2, and there is one degree of freedom for (si, S 2 , S 3 , S 4 ). □ 


The restrictions 0 and ( [Io| ) were noticed in [21] in the context of gluing tensor 
product rectangular patches with all bi(0) = —1. The restrictions are then simply 

o'i(O) = Os(0). oi(0) = <(0). 


2.4. Topological restrictions. A guiding principle for the construction of geometric 
continuous functions is that properties are equivalent to properties in the plane 
after an adequate reparameterization of the problem. Gluing two faces along an edge is 
transformed locally via such reparameterization maps, into gluing two half-planes along 
a line. Each half-plane is in correspondence with the half-plane determined by one of 
the faces and the shared edge. A natural gluing is to have the half planes on each side 
of the line. In this case, the points of one face are mapped by the reparameterizations 
on one side of the line and the points of the other face on the other side of the line. This 
implies that the transition maps keep locally the points of a face on the same side of the 
edge and thus it should have a positive Jacobian at each point of the edge. Therefore 
the hrst topological restriction that we ask for each edge r, using the canonical form 
(j^, is the following: 

Vu G [0, 1], br(M) < 0. 

When the function br is positive on the edge; the transition map identihes interiors 
of the polygons. It corresponds to two patches of surfaces virtually pasted at a sharp 
edge (i.e., at angle 0 ) rather than in a proper continuously smooth manner (i.e., at the 
angle vr). In some CAGD applications, it may be useful to model surfaces with sharp 
wing-like edges by the continuity restrictions with br > 0. But typical G^ continuity 
applications should require b^ < 0 on the whole edge to prevent this degeneration. The 
regularity property across edges, considered irrespective of orientation, is called weak 
geometric continuity [TJ §6.7] and the restriction we consider corresponds to coherently 
oriented parametrizations in [7] § 6 . 8 ]. 

Similarly, gluing the faces around a vertex 7 should be equivalent to gluing sectors 
around a point in the plane, via the reparameterization maps. Such sectors should form 
a fan around the parameter point, which can be identihed with the local neighborhood 
of the vertex 7 on the surface. Thus these sectors should not overlap. If this fan is 
dehned by vectors ui,... ,Ui 7 ’/ G (uj+i is supposed to be outside the union of the 
sectors dehned by two consecutive vectors Uj_i, for 2 < j < z), we easily check that 
the coefficients aj( 0 ), bi( 0 ) of the transition map ([^ across the edge r* at 7 are such 
that: 


( 13 ) 


Ui_i = ai(0)Uj + bi(0) Uj+i 
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or equivalently 


Ui, u 


i-i - 


0 

b*(0) 


1 

0 ^( 0 ) 


u 


*+i 


Ui 


(see also the construction in the next section 2.5). If the sector angles are less than 
71 (i.e. the sector Ui,Ui+i coincides with the cone generated by Uj,Uj+i) the condition 
that the sectors form a fan and do not overlap translates as follows: the coefficients of 
the last row of 

n( ° ' 

11 l b.(0) a,(0) 

l=J ^ 

should not be both non-negative for 1 < j ^ k < F. 

A natural way to dehne transition maps at a vertex 7 which satisfy this condition is 
to choose vectors in the plane that dehne a fan, as in Figure Then the coefficients 
aj(0), bi(0) are uniquely determined from the relations (13). 

The topological constraints could be dropped in some applications, for example, 
when modeling analytical surfaces with branching points, or surfaces with sharp wing- 
like interior edges, or with winding-up boundary. In these specihc applications, the 
compatibility Condition 2T at crossing vertex might need to be extended, to allow 
winding up of 8 , 12 , etc., crossing edges, and take into account the sharp edges. Apart 
from this kind of consideration, the topological conditions do not essentially affect our 
algebraic dimension count. 

The framework that we propose is more general than previous approaches used in 
Geometric Modeling to dehne splines (see e.g. [22l §3]) since it allows to dehne 
diherentiable functions on topological surfaces such as a Mobius strip or a Klein bottle. 

Moreover, it does not rely on the construction of manifold surfaces and atlas, but 
only on compatible transition maps. 


2 . 5 . Example. A simple way to dehne transition maps is to use a symmetric gluing 
as proposed in [T^l § 8 . 2 ] for rectangular patches. If r = (70,71) is the shared edge 
between Ui and cr2, the transition map can be of the form: 

~ (^u + 2v{'i)o{u) cos ^ - r)i(M) cos ) 

where no (resp. ni) is the number of edges at the vertex 70 (resp. 71 ). Additionally, 
if 7 o corresponds to n = 0 and 71 to n = 1, the functions a and b interpolate 0 and 
1: f)o(O) = 1, tio(l) = 0, ^i(O) = 0, ^i(l) = 1 and their derivatives of order 1 should 
vanish at 0, 1. It corresponds to a symmetric gluing, where the angle of two consecutive 
edges at 7 * is If f)o(M) and f)i(n) are polynomial functions, their degree must be 
at least 3. If f)o('u) and f)i(n) are rational functions with the same denominator, the 
maximal degree of the numerators and denominator must be at least 2 . As we will 
see the dimension of the spline space decreases when the degree increases. Thus it is 
important to construct transition maps with low degree numerators and denominators. 
See e.g. [201 El] for low degree constructions, which depend on the structure of Ai. 

A general construction of gluing data which satishes the compatibility conditions is 
as follows. 

(z) For all the vertices 7 G Ado and for all the edges Ti, ... ,tf of Aii that contain 
7 , choose vectors ui,..., G such that the cones generated by Uj, Uj+i 
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form a fan in and snch that the nnion of these cones is when 7 is an 
interior vertex. 

Compnte the transition map at 7 = (0, 0) on the edge 


i_i(0,0) = S O [ui,Uj+i] O [ui_i,Uj] O s = 


0 b^XO) 

1 a,,( 0 ) 


where S = 


0 1 
1 0 


[uj, Uj+i] is the matrix which colnmns are the vectors 


u,, u 


■i’) *^2+1 5 


\Ui,Uj 


is the determinant of the vectors Uj, u. and 


ar,( 0 ) = 


1 ^ 2 — 1 : I 


■ 6 ^,( 0 ) = 


1 ^ 2 , U2-I-I I 


(u) For all the edges t E A4i, dehne the rational fnnctions = ^, b^ = ^ on the 
edges r by interpolation as follows: if there is no crossing edge in then 
a linear interpolation of the valne at the vertices is snfficient. If 71 ,..., 7 ^ 
is a seqnence of crossing vertices, and 7 o 7 i, 7 i 72 , • • • ,7n7n+i is a seqnence of 
edges passing “straight” throngh them, we can choose linear glning data on 
one edge, and qnadratic data on the remaining edges of the seqnence so that 
the constraints (|^ and ( 10 ) are satished. 

Therefore, for general meshes, glning data which satisfy the compatibility condition 
and the topological condition can be constrncted in degree ^ 2 . 


3. Spline space on a topological surface 


The main object of onr stndy is the space of fnnctions on the topological snrface 
Al, which are differentiable and piecewise polynomial. Such functions are called spline 
functions on A4. Let 7^(cr) = be the ring of polynomials in the variables 

{ua, Va) attached to the face a. A spline function / is dehned by assigning to each face 
<J G M .2 a polynomial fa G 7Z(a), and by imposing the regularity conditions across the 
shared edges. 

We also consider rational gluing data on the interior edges r G Adi: 


(15) 


cIt-(mi) 


Ct(Ui) 


and bT-(Mi) 


brjui) 

Cr(ui) 


with ar(ui), bT-lui) and Cr(ui) polynomials in the variable ui, where br(Mi) and Cr(ui) 
do not vanish on r (i.e. for Ui G [0,1]). As b(ui), c(ui) do not vanish on r, the transition 
map 00 - 2 , 0-1 is a (F^-diffeomorphism in a neighborhood of the edge r = (ri,r 2 ) between 
(Ti and (Ti. The polynomial vector [a,-, c.^] is also called hereafter the gluing data of 

the edge r. We assume hereafter that the transition maps satisfy Conditions 2.6. 2.1 
and all crossing vertices of Ai have 4 edges. 

We can now dehne the space of splines on Ad: 


Definition 3.1. We denote by 5^(Ad) the (M-linear) space of differentiable functions 
on the topological surface Ad which are dehned by assigning polynomials to the faces 
a E A4 2 satisfying the constrains ([^. More precisely, 

S\M) := {/ G ©aeA^27^(c^) | Jy(fai) = d.^(/o 2 o V 7 G r = cti H da}, 

where Jj is the jet or Taylor expansion of order 1. 
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A spline / G gives a piecewise-polynomial map, defined on every face of A4, 

and the jets of order 1 coincide on the shared edges. This definition implies that a 
spline function / G (Ad) is on a neighborhood of a shared edge r = cxi fl <72 if we 
use the re-parametrization o (()a 2 ,( 7 i- 


Dehnition |3.1| can be directly extended to splines of any order r but in this paper 
we only consider r = 1. 


3.1. Polynomials on faces. On each face a G Ad 2 , we consider polynomials of degree 
bounded by fc G N. 

If cr G Ad 2 is a triangle (P, Q, R), we denote by the finite dimensional vector 

space of polynomials in with total degree bounded by k. 

After a change of coordinates, we may assume that the coordinate function u = 
satisfies Ua{PR) = 0 and Ua{Q) = 1, while v = satisfies Va{PQ) = 0 and Va{R) = 1. 
Introducing tc = 1 — m — n, we can express any polynomial in TZkicr) as a homogeneous 
polynomial of degree ^ k in the barycentric coordinates u,v,w using the Bernstein- 
Bezier basis: 



k\ 


i\j\{k - i-j)\ 


u 




for 0 ^ i ^ fc. We verify directly that for a function 


/(m, v)= ^ Cijbfj{u, v,l-u-v) 

O^i+j^k 


expressed in this basis, we have 


(16) /(0,0) = co,o, 

duf{0, 0) = /c(ci,o - Co,o), dyf{0, 0) = /c(co,i - Co,o), 

dudyf{0, 0) = k{k-l) (ci,i - ci,o - co,i + co,o)- 


If a G Ai 2 is a rectangle {P,Q,R,S), we will denote by lZk{,cr) the hnite dimensional 
vector space of polynomials in M[Mo-,fcr] with partial degree in and bounded by 
k, where u = is chosen such that Ua{PS) = 0, Ua{QR) = 1, and n = Vo- is chosen 
such that v„{PQ) = 0, Va{RS) = 1. Introducing h = 1 — h = 1 — n, we can express 
any polynomial function of TZk{o') as a bi-homogeneous polynomial of degree k inu,u 
and degree k in v, h, using the tensor product Bernstein-Bezier basis 


h-/u,u,v,v) = — 


k\k\ 


i\j\{k - i)\{k - j)\ 


yi yk-i yjyk-j_ 


for 0 ^ i ^ fc, 0 ^ ^ fc. We verify directly that for a function / = j^k 

expressed in this basis, we have 


(17) /(0,0) = co,o, 

duf{0, 0) = fc(ci,o - co,o), dyf{0, 0) = k{co,i - co,o), 
dudyf{0,0) = (ci^i — Ci^o ~ Co,i + Co,o)- 

The finite dimensional vector space of spline functions / = ifa)( 7 eM 2 ^ ‘5(A1) of degree 
bounded by fc G N on each face (/o- G 7lk{o')) and of regularity r is denoted SJl(M) or 
simply Sk(M) when r = 1. 
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3.2. Taylor maps. An important tool that we are going to nse intensively is the 
Taylor map associated to a vertex or to an edge of A4. 

Let 7 G Mo be a vertex on a face a G M 2 belonging to two edges r, r' G Mi of 
a. We dehne the ring of j on a by 7 ^°" ( 7 ) = TZ{a) f {il, i^,) where is the ideal 

generated by the squares of ir and ir', the equations ir{u,v) = 0 and ir'{u,v) = 0 are 
respectively the equations of r and r' in 7^(a) = M[M,n]. 

The Taylor expansion at ■y on a is the map 

T; : / G 7^(a) ^ / mod in 7^"(7). 

Choosing an adapted basis of one can dehne by 

Tpf) = [/(7). a„/(7), 9./(7), 9A/(7)]. 

The map can also be dehned in another basis of in terms of the Bernstein 

coefficients by 

7 (/) = [C(,,„(/), Ci,„(/), c„,,(/), c,,i(/)] 

where co,o, ci,o, co,i, ci^ are the hrst Bernstein coefficients associated to 7 = (0, 0). 

We dehne the Taylor map on all the faces a that contain 7 , 

T-y:f = (/.) G ©.7^(a) ^ G ©.^^7^"(7). 

Similarly, we dehne Tq as the Taylor map at all the vertices on all the faces of M. 

For an edge r G Adi on a face a G M 2 , we dehne the ring of t on a by VTir) = 
7?.(ct)/(£^) where ir{u,v) = 0 is the equation of r in 7l{cr) = ]R[M,n]. The Taylor 
expansion along r on a is dehned by 

T; : / G 7^(a) ^ / mod (£?) in 
and the Taylor map on all the faces a that contain r is given by 

Tr-.f = (/.) e ©.7^(a) ^ {TfiU)) G ©.^.7^"(7). 

Similarly, we dehne Ti as the Taylor map along all the edges on all the faces of M. 

4 . SPLINES ALONG AN EDGE 


To analyze the constraints imposed by gluing data along an edge, we consider hrst 
a simple topological surface J\f composed of two faces cxi, (J 2 glued along an edge r. 

A spline function / G on Af is represented by a pair of polynomials / = (/i, / 2 ) 

with fi G Tl{ai) = M[ni, Vi] for z = 1, 2. 

By a change of coordinates, we assume that the edge r is dehned by ui = 0 and 
Ml G [0,1] in (Ti and by M 2 = 0 and M 2 G [0,1] in a 2 . 


4.1. Splines and syzygies. With the transition map 00-2,0-1 dehned by the rational 
functions a = ^ and b = — as in (15), the diherentiability Condition (4) along the 
interior edge r becomes 


a(Mi)A(Mi) + 6(Mi)i?(Mi) + c(mi)C(mi) = 0, 


where 

A(mi) = ^(0,Mi), 5 (mi) = ^(0,Mi), (©(mi) =-^(mi,0). 

OV 2 OU 2 OVi 

Thus, the (©^-smoothness condition along an interior edge is equivalent to the condition 
on (A, B, C) of being a syzygy of the polynomials a(Mi), 6 (mi), c(mi). 
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More precisely, a spline (/i,/ 2 ) on M is constructed from a syzygy {A,B,C) of 
a, b, c by defining: 


(18) 

(19) 


fUi 


fi{ui,vi) = co+ / A{t)dt - viC{ui) + v^Ei{ui,vi), 


rV2 


f2{u2, V2) = Co + / A{t)dt + U2B{v2) + M2-E2(m2, V2), 


where cq G M is any constant, and Ei, E 2 are (any) polynomials in Vj\ for i = 1,2, 
respectively. 

We will use this representation for the splines on W to compute the dimension of 
the space of splines SKAf), see Proposition 4.6 below. Before, we introduce some 
notation, both for the proof and the dimension formula. 


The module of syzygies of a{ui),b{ui),c{ui) over the ring M[mi] is denoted by Z = 
Syz(a, b, c). For (^4, B, G) G Z, the maximum of the degrees, max(deg A, deg B, deg G) 
is called the coefficient degree of the syzygy. 

Each of the faces ci and (J 2 in Af can be a triangle or a rectangle. Let us denote by 
En the number of rectangles and by Fa the number of triangles in Af. 


Definition 4.1. As before, let Ui, a 2 be the faces of Af. We dehne 

m = min (Fa (ai). Fa ( 0 - 2 )), 

where Effiaf) = 1 if cij is a triangle and 0 otherwise. For the polynomials a,b,c E M[mi] 
dehning the gluing data along the edge r, let n = max(deg(a), deg(6), deg(c)), 

da = n + l, dh = n + Ea{<J 2 ), and dc = n + Effiai), 


and 



if min((ia — deg(a), db — deg(6), dc — deg(c)) = 0 and 
otherwise. 


By the formulas (18) and (19) representing a spline (/i,/ 2 ) G SKAf), let us notice 
that we need to consider syzygies {A, B, G) of a, b, c such that deg(A) ^ k — 1, deg(F) ^ 
k — Fa(( 72), and deg(F) ^ k — Effiai). The reason is that /* is of bidegree at most 
(fc, k) if CTj is a rectangle and of total degree at most k if cxj is triangle. 


Definition 4.2. For fc ^ 0, we will denote by the vector subspace of Z of syzygies 
of {a,b,c) dehned as the set 

Zk = {{A, B,G) e Z: deg(y4) ^k - 1, deg(F) ^ k - Fa(< 72 ), 

and deg(F) ^ k — Fa((7i)}. 

Let us consider the map 
(20) Qr-.Z^ S\Af) 

( pui rv2 \ 

J A{t)dt — ViG{ui),j A{t)dt + U 2 B{v 2 )\ ■ 

By construction, we have 0T-(Zfc) C Sl{Af). 

The dimension of Z^, as a vector space over M, will be deduced from classical results 
on graded modules over S = M[mo,Mi]. We will study the module Syz(a,6, c), where 
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a,b,c E S are the homogenization of a, b and c in degree da, df, and d^ respectively. The 
elements in Syz(d, b, c) in degree n + k will precisely lead to the syzygies in Z^- 

Lemma 4.3. For polynomials a,b,c E with b,c 0, gcd(a, 6, c) = 1 and Z = 

Syz(a, b, c) as defined above, 

(i) Z is a free Mfufi-module of rank 2. 

(m) The module Z is generated by vectors {A 2 , 82 , 02 ) of coefficient 

degree /i and u = n — jj, + 1 + F/^ — e — 2 m where fi is the smallest possible 
coefficient degree. 

(m) For k eN, the dimension of Z^ as vector space over M is given by 
dim Zk = {k — pi — m + 1)+ + {k — n + pi + m — Fa + e)+ 
where = max(0, t) for t eTj. 

{iv) The generators {Ai,Bi,Ci), {A 2 , 82 , 02 ) of Z can be chosen so that 
{a, b, c) = {B 1 O 2 — B 2 O 1 , O 1 A 2 — O 2 A 1 , A 1 B 2 — A 2 B 1 ). 


Proof. We study the syzygy module Z = Syz(a, b, c) using results on graded resolutions. 
For this purpose, we homogenize a, b and c in degree da = n + 1, dh = n + Fa (< 72 ), and 
dc = n + Fy\{ai), respectively, where F/^{ai) is as in Dehnition 4.1 Let Uo,Ui be the 
homogeneous coordinates, and d, b, c the corresponding homogenizations of a, b, and c. 
We consider the module of homogeneous syzygies Syz(d, b, c) over the polynomial ring 
S = ]R[uo, ui]. 


Claim 4.4. For any k ^ 0, the elements in Zk are exactly the syzygies of degree n + k 
in Syz(d,6, c) after dehomogenization by setting uq = 1. 


Proof. It is clear that if Ad + Bb + Oc = 0, then by dehomogenization taking mq = 1, we 
get a syzygy {A, B, O) of (a, b, c). Moreover, if deg(yla) = deg{Bb) = deg(Cc) = n + k, 
then deg(A) = k — 1, deg(F) = k — F/^{a 2 ) and deg(C) = k — Fa(cti). It follows that 

{A,B,0)EZk. 

On the other hand, any syzygy {A,B,0) E Zk is given by polynomials that satisfy 
the conditions in Dehnition 4.2 Thus maxjdegA, deg F, degC} ^ k, and since n = 


max{dega, deg 6, degc} then we may consider the homogenization of the polynomial 
Aa + Bb + Oc in degree n + k. These polynomials satisfy 


0 = u^+'^iAa + Bb + Oc){ui/uo) 

= Aa{ui/uo) + Uq Bb[ui/uo) 

It is easy to check that 

A = Uo~^A{ui/uo), B = u4^^^‘'''^B{ui/uo), C = 

are all polynomials in M[mi,mo]) and dehne a syzygy of d,b,c of degree n + k. Let us 
also notice that the polynomials 

a = Mo’'"^a(ui/uo), b = u4^^^'^^'^b{ui/uo), and c = 


are precisely the homogenization of a, b, c in degree da, db, dc, respectively. □ 
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As gcd{a, b, c) = 1, we have gcd{d, b, c) = Uq if e = 1, and gcd{d, b,c) = 1 otherwise. 
Let / = (a, 6, c) be the ideal generated by d,b,c in S. If gcd{d,b,c) = 1 then there 
exists to E N such that Wt ^ to, It = ( mo , mi )* and in that case, dim^{S/I)t = 0 for 
t sufficiently large. It follows that the Hilbert polynomial HPs/i of S/I is the zero 
polynomial. 

For the second case, namely if gcd{a,b,c) = uo, since gcd(a, 6, c) = 1 then the 
polynomials d/uo, b/uo and c/uo have gcd equal to 1. Hence there exists to ^ N such 
that Vt ^ to, It = Uo{uo,Uiy~^. In this case dim]^{S/I)t = 1 for t sufficiently large, 
and it follows that the Hilbert polynomial HPs/i is the constant polynomial equal to 

1 . 

Then the exact sequence 

0 ^ I ^ S ^ S/I ^0 

implies that 

(21) HPi{t) = HPs{t) - HPs/i{t) = - e, 

where HPm is the Hilbert polynomial of the module M. 

By the Graded Hilbert Syzygy Theorem, we get a resolution of the form 

0 S{-d^) © • • • © S{-dL) ^ S{-da) © S{-db) © S{-d,) ^ 0 . 


Notice that this resolution is not necessarily minimal. Since this is an exact sequence, 
then the Hilbert polynomial of the middle term is the sum of the other two Hilbert 
polynomials, and applying (21) we get 

3t — (dfl + dfo + dc) + 3 = (t — di + 1) + ■ ■ ■ + (t — d/, + 1) + (t + 1) — e. 


It follows that L = 2 which proves (i). Furthermore, we have that the degrees di and 
d 2 of the syzygies satisfy di + d 2 = d^ + df, + dc — e. 

The matrix A representing A is a 3 x 2 matrix 

/A, A2 \ 

Ml ^2 

V C2 J 


the hrst column corresponding to the generator of degree di and the second of degree 
d 2 . These two syzygies correspond to vectors of polynomial coefficients of degree /i = 
di — min(da, df,, dc) and z/ = d 2 — m.m{da, db, dc)■ By Definition 4.1, min(da, d;,, dc) = 
n + min(l, Fa(cti), Fa(ct 2 )) = n + m, and also d^ + d;, + dc = 3n + Fa + 1. Let us 
assume that di ^ d 2 , then ft is the smallest degree of the coefficient vector of a syzygy 
of (d, b, c), and z/ = n — /i + I + Fa — e — 2 m. 

By exactness, the two columns of A generate Syz(d, b, c). The dehomogenization (by 
setting Uo = 1) of the syzygies in Syz(d,6, c) leads to syzygies of {a,b,c) over M[mi]. 
In particular, it is straightforward to show that the dehomogenization [At, Bt, C/) of 
{Ai,Bi,Ci) for i = 1,2 generate Z = Syz(a, 6, c) as a module over M[mi]. This proves 
(^^)- _ 

the space is in correspondence with the space of homogeneous 


By Claim 4.4 


syzygies of degree n+k, which is spanned by the multiples of degree n+k of {Ai,Bi, Ci) 
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and (^ 2 ,-02, ^' 2 )- Therefore, 

Aim Zk = {n + k - di + l)+ + {n + k - d2 + l)+ 
= {k — ^ — m + 1)+ + {k — v — m + 1)+ 


with i/ = n — /i + l + Fa — e — 2m. This proves (Hi). 

The point (iv) is a conseqnence of Hilbert-Bnrch theorem. More details on this proof 
can be fonnd in [6l Chapter 6, §4.17]. □ 


Definition 4.5. For an interior edge r in the topological snrface Ai shared by the faces 
(Ti, (72 with glning data [a,-, b^, c,-], we denote by Mr the topological snrface formed by 
the cells ui, a 2 glned along the edge r with the same glning data. Let fir be the 
smallest coefficient degree among the two generators of the modnle Z = Syz(aT-, br, Cr). 
Let Ur = Ur — fir + I + “ 6,- — 2771,- denote the complementary degree, where 

Ur = max(deg(ar), deg(cr), deg(cr)), mr = mm{FA{ai), Fa{o- 2 )), ^r = min(?7,^ + 1- 
deg(a^),n^ + FA(cr 2 )-deg(6^),n^ + FA(cri)-deg(c^)) FA(r) = FA(ai) + FA(a 2 ). The 
corresponding basis of the syzygy modnle Z of [a,-, br, Cr] is called the pT-basis. 

This constrnction allows ns to determine the dimension of S\{Mr). 

Proposition 4.6. For F^ (resp. Fa) the number of rectangles (resp. triangles) of Mr, 

dimSk{Mr) = 1 + - l)Fn + - k)FA + dr{k) 

where dr{k) = {k — fir — rrir + 1)+ + {k — Ur + fir + rrir — Fa + er)+ 


Proof. Since the only constraints satished by the spline fnnctions in S{Mr) are the 
glning conditions along the edge r, the nnmber of linearly independent splines on Mr 
can be easily connted by nsing (18) and (19), and the linearly independent terms in the 
Bernstein-Bezier representation of the polynomials fi and /2 that conform a spline. 

The glning data and the smoothness along the edge r impose conditions on the terms 
in fi and /2 which are linear in vi and V 2 , respectively. Thns, the dimension of the 
space of splines on Mr of degree exactly 1 in vi and V 2 , is given by dimZ^ = dr{k). 
The formula for dr follows from Lemma 4.3, by considering Z = Syz(ar, &t, Cr), where 
Or, &r, Cr dehue the gluing data along r. □ 


4.2. Separation of vertices. We analyze now the separability of the spline functions 
on an edge, that is when the Taylor map at the vertices separate the spline functions. 

Let / = (/i, / 2 ) G 77.((Ti) © 77.((T2) of the form /*(«*, Vi) =Pi + qiUi + qi Vi + Si UiVi + 
riuf F fivf F ■ ■ ■. Then 

MU) = [PD9D?lAl,P2,g2,g2,S2]. 

If / = (/i) f 2 ) G Sl{Mr), then taking the Taylor expansion of the gluing condition (|^ 
centered at = 0 yields 

(22) ^i + SiMi = (ci(0) + Cl ( 0 ) 7 X 1 + ■ ■ ■ ) (^2 + 2 72 Ml + ■ ■ ■) 

+ (b(0) + b^(0)Mi + ■ ■ ■) (q'2 + S 2 + ■ ■ ■) 
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Combining (22) with ([^ yields 


Pi = P2 

qi = q 2 


ri = r 2 

qi = a(0)g2 + l’(0)g2 

Si = 2 ci(0) r2 + b(0) S2 + ci^(O) q2 + b^(0) q2- 


Let 7 /( 7 ) be the linear space spanned by the vectors [pi, gi, gi, Si,p 2 , g 2 , ? 2 , S 2 ], which 
are solution of these equations. 

If a(0) 7 ^ 0, it is a space of dimension 5 otherwise its dimension is 4. Thus dim 7 /( 7 ) = 
5 - c^( 7 ). 


Proposition 4.7. For k ^ + m-r + 1, (AC)) = 'HH). Its dimension is 

dimT,(5,HAC)) = 5 -c,(7). 

Proof. Let ^(7) = T.y( 5 ^(A/"r))- By construction (^(7) C ^(7). We are going to prove 
that for k ^ Ut- + rUr + 1, G{'j) and ^( 7 ) have the same dimension and thus are equal. 
By the decompositions (18) and (19), the elements of T^{Sl{Mr)) are of the form 

[co, A(0), -C'(O), -C"(0), Co, 5(0), A(0), B'i 
and [A, 5,(7] G Z^. By Lemma 4.3 


where Cq G 


an element of is of the form 
[A, B,C] = P [Ai, Bi, Cl] + Q [A 2 , B 2 , C 2 ] with P,Qe M[m], deg(P) ^ fc — p,- — 771 ^- and 
deg((5) ^ k — Ut- — nir- By removing the repeated columns, reordering and changing 
some signs, we see that ^( 7 ) = T.y(5^(AC)) is isomorphic to the space spanned by the 
elements 


(23) 


flip) 


0 

0 

0 

0 


r 1 

dujiil) 

^^ 2 / 2 ( 7 ) 


0 Ai(0) ^ 2 ( 0 ) 0 0 

0 Pi(0) Psjo) 0 0 


Co 

P(0) 

Q(o) 

P'(0) 

Q'(o) 

-dvjiil) 

du2dvj2{l) 


0 Ci(0) (© 2 ( 0 ) 0 0 

0 B[{0) P'(0) Pi(0) P2(0) 


. -duAJi{-r) . 


L 0 C((0) C'(0) CCO) ^2(0) J 



for P,Qe M[m] with deg(P) ^ k — fir — nir and deg(Q) ^ k — Hr — nir. Let us assume 
that k ^ Ur + nir + 1 so that k — fir — nir ^ k — Vr — nir ^ 1- 

As A 1 P 2 — A 2 P 1 = c and Ai( 0 )P 2 ( 0 ) — A 2 ( 0 )Pi( 0 ) = c( 0 ) C 0) we deduce that 
[Pi(0),P 2 ( 0 )] C [0)0] cind that dimG( 7 ) ^ 4. 

If Ct-( 7 ) = 0, then a(0) = Pi( 0 )C 2 ( 0 ) — P 2 ( 0 )Ci( 0 ) C 0 and dimG( 7 ) = 5 = 
5 — c.,-( 7 ) = dim 7 ^( 7 ). 

If Ct-( 7 ) = 1, then a(0) = Pi( 0 )C 2 ( 0 ) — P 2 ( 0 )Ci( 0 ) = 0 and dimG( 7 ) = 4 = 
5 — c.,-( 7 ) = dim 7 ^( 7 ). 

In both cases, we have dimG( 7 ) = dim 7 /( 7 ), which implies that ^( 7 ) = ^( 7 ). This 
completes the proof of the proposition. □ 


If 7 ' is the other end point of r, we have a Taylor map for each 7 and 7 ', that we 
join together. Let 

(24) T.,,y ;P(ai)©p(a 2 ) ^ 

/=(/l,/2) ^ (T,(/),Ty(/)) 

and let G{t) = (5^(AC))- 
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Proposition 4 . 8 . For k ^ Ut- + rriT- + 4 :, we have (iS^( 7 V"t-)) = 'H{'y')) and 

dimr^,y(5fc(7V;)) = 10 - c^(7) - c^(7'). 

Proof. By a change of coordinates, we can assnme that the coordinates of 7 (resp. 7') 
are (0,0) (resp. (1,0)) in ui and (0,0) (resp. (0,1)) in (72. 

Similarly to the proof of the previons proposition, Ty^SKAfr)) is spanned by the 
vectors 

[co+ f A{u)du,A{l),-C{l),C\l),co+ [' A{u)du,B{l),A{l),B\l)] 

Jo Jo 

for Co G M and [A, B,C] = P [Ai, Bi, Ci] + Q [^2, B2, C2] G with deg(P) ^ k — fXr — 
rrir and deg(Q) ^ k — Ur — 

For k ^ Ut + rUj. + 4, we can find polynomials P = po(l ~ + 2vJ‘) +pi{u — 2vJ — 

u^) + P 2 V?{1 — uY, Q = qQ{l — ?>v? + 2uY + qi{u — 2v? — v?) + q 2 U^Y- — uY, of degree ^ 4 
snch that P(0) = po, -P'(O) = pi, Q(0) = go, Q'(0) = gi, P(l) = 0, P'(l) = 0, Q(l) = 0, 
Q'(l) = 0 and f^(PAi + QA 2 ){u)du = —Cq. 

This implies that ('H(7),0) C G{t). By symmetry, we also have (0,^(7')) C G{t). 
By constrnction G{t) C therefore we have G{t) = (^(7),^(7')) and 

dimG(r) = dim'P(7) + dim7/(7'). We dednce the dimension formnla from Proposition 

SjI □ 


Definition 4.9. The separability s(r) of the edge r is the minimal k snch that 

T.ysKMV) = {T,{si{K)),Tysi{K))). 

Remark 4.10. The bonnd + m,- + 4 ^ s(r) is not necessarily the minimal degree 
of separability. Separability can be attained as soon as dr{k) ^ 9 — c.,-( 7 ) — Ct-( 7 '). 

4.3. Decompositions and dimension. Let r G Ali be an interior edge r shared 
by the cells (7i,(72 G M. 2 . Let Ki = {vf) fl Pfc((7i) and K 2 = (n^) F T^k{a' 2 ) be the 
polynomials of Pfc(cri) (resp. Pfc((72)) divisible by vf (resp. ul). Let L be the snbspace 
of polynomials of Pfc(o'i) ©”^^(^ 2 ) spanned by the Bernstein basis fnnctions on ai and 
(72, which are not divisible by vf or and let til be the projection of 7lk{o'i) ©Pfc(o' 2 ) 
on L along (Pi,0) © (0,P2)- The fnnctions in L are said to have their support along 
T. By constrnction, we have Pfc((7i) © Tlk{o' 2 ) = (Pi,0) © ( 0 ,^ 2 ) © L. The elements 
of (Pi, K 2 ) are obvionsly in SI{Nt) since they vanish at the order 1 along r. 

Let hFfc(r) = niiPiYZk)) where 0,- is defined in (20). Notice that W4(r) C Sl{A/r) 
since kerTr^ C Sl{J\fr). Moreover, since kerTr^ does not intersect ©^(Z^) and 0.^ is 
injective, the spaces WYt), 0T-(Zfc) and have the same dimension. Therefore, we 


have dim(lFfc(r)) = dr{k) and Wkir) Y {0} when k ^ + m (Lemma |4.3| (hi)). 

From the relations 

(25) 


and ( |19[ ), we dednce the following decomposition: 
Sl{K) = (^ 1 , 0) © (0, P 2 ) © Mu © IFfc(r) 


where u = 7rL((l,l)). The snm of these spaces is direct, since the snpports of the 
fnnctions of each space do not intersect. 

The map defined in (24) indnces the exact seqnence 

0 ^ /Cfc(r) ^ Sl{Nr) ^ G{t) ^ 0 

where /Cfc(r) = kerPy^ G{t) = T.y^y{Sl{A/r))- If is clear that {Ki,K 2 ) C ICkir). 
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Definition 4.11. For an interior edge r G M°, let = ker(T^_y)nlFfc('r) be the set 
of splines in Sl{J\fr) with their support along r and with vanishing Taylor expansions 
at 7 and 7'. For a boundary edge r' = (7,7'), which belongs to a face a, we also define 
Tfc(r') as the set of elements of 7lk{o') with their support along t' and with vanishing 
Taylor expansions at 7 and 7'. 

Notice that the elements of Tfe(r) have their support along r and their Taylor expan¬ 
sion at 7 and 7' vanish. Therefore, their Taylor expansion along all (boundary) edges 
of A/'t distinct from r also vanish. 


Lemma 4.12. For an interior edge r E we have JCkir) = (iFi, 0)©(0, iF 2 )©i/fc('r)- 

Proof. As (A'ljO), ( 0 , 1 ^ 2 ) C kerT.y^Y = and JCkir) fl (IFfc(r) ©Mu) = ICkir) fl 

Wkir) = Tfc(r), we have 

/Cfc(r) = (i^i,0)©(0,iF2)©((lFfc(r)©Mu)n/Cfc(r)) 

= (Ai, 0) © (0, 1 ^ 2 ) © ^^fc('7‘). □ 


Corollary 4.13. For an interior edge r G A4° and for k ^ s(r), the dimension of 
£k{r) is 

dimSkir) = d^ik) - 9 + Ct(7) + Ail')- 


Proof. By Lemma 4.12, we have 

dimTfc(r) = dim/Cfc(r) — dim Ai — dim A 2 . 
As ICkir) is the kernel of and G(r) is its image, we have 

dim/Cfc(r) = dim SliPfr) — dimG(r). 


As dim(lFA;(r)) = drik), we deduce from the decomposition (25) that dim (A/'t) = 
1 + dr{k) + dimAi + climA 2 . Using Proposition 4^ G'(r) = {'H{i),'H{i')) and we 
obtain 


dim£/fc(r) = dimS'^(A/'T) — dimG(r) — dim Ai — dim A 2 
= drik) - 9 + Ct(7) + Ct(7')- 


□ 


Remark 4.14. When r is a boundary edge, which belongs to the face cxi G AI2, 
we have Ski^fr) = T^kio'i), Akir) = Ai © Tfc(r) and for k ^ 5(r) = 3 + FA(ai), 
dimG(r) = 8 and dimTfc(r) = k + l + ik + l — F^iai)) — 8 = 2 k — F/^iai) — 6. 

Notice that this is also what we obtain if we attach a virtually rectangular face along 
T with constant gluing data: n = /i = 0, m = 0, e = 2, Ct(7) = Ct(70 = 0 

drik) = 2k + 3 - F^iai), 
so that dimTfe(r) = drik) — 9 + Ail) + Ail')- 


5. SPLINES AROUND A VERTEX 

We consider now a topological surface O composed of faces cxi,..., <7^ G O 2 sharing a 
single vertex 7, and such that cXi and cxi+i share the edge Tj+i = (7, (5j+i). In particular 
Tj, Tj+i are the two edges of Uj containing the vertex 7. The number of edges containing 
7 is denoted F'. All the vertices of O different from 7 are boundary vertices. The vertex 
7 is an interior vertex, iff ap and Ui share the edge Ti. In this case, we identify the 
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indices modulo F and we have F' = F, otherwise we have F' = F + 1. The gluing 
data for the interior edge r* is a* = —, b* = —. 

The coordinates in the ring 7l{ai) are chosen so that the coordinates of 7 are (0,0) 
and Ti is dehned by Vi = 0, Ui E [0,1] and by Mj_i = 0, G [0,1] in 7?.((Tj_i). The 
canonical form of the transition map at 7 across the edge r* is then 


: {u, v) 


Vi bi{ui) 

Ui T Vi Cli {Ui 


Let / = (/i)i=i,,..,F € The gluing condition (|^ implies that the Taylor expansion 

of fi at 7 is of the form 


fi{ui, Vi) = p + Qi Ui + qi+i Vi + Si UiVi + Ti uj + Tj+i H- 

for p,qi,Si,ri eR, i = 1,... ,F (see Fig. |^. 



Figure 3. Taylor coefficients around a vertex. 

By a computation similar to (22), Condition (|^ implies that 

(26) qi+i = 0^(0) qi + bi(0) qi_i {i = 2,...,F) 

(27) Si = 2 aj(0) r* + bi(0) Sj_i + a((0) (ji + b((0) (ji-i {i = 2,...,F) 

Let 77 ( 7 ) be the vector space spanned by the vectors h = \p,qi,... ,qFi, si,..., sp] 
for h' = [p,qi,... ,qpr,si,... ,SF,ri,... ,rFi] a solution of the linear system (26), (27). 

Proposition 5.1. 

dimH( 7 ) = 3 + F( 7 ) - ^ c^( 7 ) + £+( 7 ) 

r97 

where F = F{^) is the number of faces around the vertex 7. 

Proof. Notice that 77 ( 7 ) is isomorphic to the projection of the solution set of system 


(26), (27) on the space of the variables [p, q, s] = [p, gi, ..., gi?/, si, ..., s^]. 

The solutions in q = (gi,... ,qF') of the first set of equations satisfy the induction 
relations 

(li \ = f 0 1 W 7*-i 

gi+i J I bi(0) ai(0) J \ qi 


for i = 2 ,..., F. 


As we have the compatibility condition 2.6 at an interior vertex, the solutions of (26) 


span a linear space of dimension 2 , parametrized for instance by gi, g 2 . 


The system (27) is formed by linearly independent equations which involve and 


g,, Sj when 0 ^( 0 ) 7 ^ 0 and by equations which only involve Sj, Si_i and qj when ai(0) = 0. 

Therefore the projection of the solution set of (27) on the space corresponding to 
the variables [p, q, s] is defined by the equations which only involve s*, Sj_i and g*, gj_i 
when aj( 0 ) = 0 . 
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If one of the edges around 7 is not a crossing edge, then the codimension of this 
space is XItb-t with the convention that Ct-( 7) = 0 if r is a boundary edge. 


If all the edges around 7 are crossing edges, then the compatibility conditions 2.7 


at 7 imply that one of these equations is dependent from the other. Therefore, the 
codimension of this space is ~ ^+(7). 


By intersecting it with the solution space of ( 26 ), we deduce that the dimension of 
77(7) is precisely given by 


1 + 2 + F 


t97 


-(7) + c+(7). 


□ 


Let = 11(737 ^7 Taylor map at 7 on (9 and let Tqq = be the 

Taylor map along all the boundary edges which do not contain 7. 

For fc G N, we dehne Vfc(7) = kerTgo fl Sl{ 0 ) the set of spline functions on O 
which vanish at the hrst order along the boundary edges (which do not contain 7). 

Proposition 5.2. For k ^ maxj=i^...^i7’(s(ri)), T..y(Vfc( 7 )) = T-LAl)- 


Proof. By construction, the elements of Vfc(7) satisfy the equations ( 26 ), ( 27 ). This 
implies that (12^(7)) cFLAl)- 

Consider an element h = (hi,..., hp) G 'HAl)- By Proposition 4.8, for k ^ s(rj), 
there exists (/*,/*) G such that T.^(/j,/j) = (h*, hi_i) and TsXfhfi) = 0. Let 

Vi = 0 (resp. Ui-i = 0 ) be the equation of r* in cTj (resp. cTj_i). As for any polynomials 
p G AA n TZk{(Ji),q G {ui_iY n T^{p, q) = 0, we can assume that (/*, /*) has 

its support in 

By construction, we have T^*(/i) = T^*(/j+i) = h*. Thus, there exist Pi G TlkAi) 
supported in VfAp) + 77 .'^‘(©+i) such that Tf^A9i) = fi, = /*+i- ^ is con¬ 

structed by taking the coefficients of fi on TlF^iji) and those of /i+i on 77 .°'*(rj+i), the 
coefficients in HV/Af+A coinciding (see Fig. |^. As Tff{fi) = Tff{gi) = 0, 

^7+i(/i) = ^^1+1 (hi) = 0 and Qi is supported in n^'{Ti)+ n^AF+i)^ we have Tf^igi) = 0 
for any edge r of the face cTj, which does not contain 7. 


Ti+l 



Figure 4 . Lifting the edge functions. 

Let g = [gi,... ,gF] G (Bai^'y'J^kio'i). By construction, g vanishes at the hrst order 
along all the boundary edges of O, which do not contain 7. Moreover, Tr^(g) = (/«, fi) G 
SKAfrJ, thus g satishes the gluing conditions along the edge r*. We also have Tr(g) = 0 
for any edge r, which does not contain 7. Thus g satishes the gluing conditions along 








22 


B. MOURRAIN, R. VIDUNAS, AND N. VILLAMIZAR 


all the edges and its image by Tqq vanishes, i.e. g G Sl{0) fl kerTgo = Vkis)). 
By constrnction, T-y(g) = h. This shows that T-Lij) C T^(Vfc( 7 )) and conclndes the 
proof. □ 


6. SPLINES ON A GENERAL MESH 

We consider now a general mesh Ai with an arbitrary nnmber of faces, possibly with 
bonndary edges. 

We denote by Tq = n 7 eAio Taylor map at all the vertices of M. and l-i = 

To{Sl{A4)). We have the following exact seqnence: 

0 ^}Ck^Sl{M)^G^0 

where = kerTo fl Sl(A4) and G = To(Sl(A4)). Let s* = max{s(r) | r G Aii). We 
have s* ^ max{z/T- + + 4 | r G Aii} 


6.1. Splines at a vertex. Let 7 G Aio be a vertex of Ai and let be the snb-mesh 
associated to the faces of Ai which contain 7 . Let Vk{l) be the set of spline fnnctions 
in Sl{Ai) snpported on the faces of O^, which vanish at the hrst order along the edges 
that do not contain 7 . 


Proposition 6.1. For k ^ 5*, ToiSliM)) = and 

dimTo(5fc^(Wf)) = ^ (^( 7 ) + 3) - ^ ^c^( 7 )+ ^ c+( 7 ), 

'y&Mo 7SAI0 ''■37 7 eAIo 


where F{^) is the number of faces of Ai that contain the vertex 7 G Aio. 

Proof. By proposition |5.2 for k A s* the image of Va;( 7 ) by Tq is (^( 7 ), and Ty (Va:( 7 )) = 
0 for any other vertex 7 ' 7 ^ 7 . 

This shows that To(iS^(Al)) = 117 ^( 7 ) = 117 dednce the dimension 
formula from Proposition |5.1[ □ 


6.2. Splines on edges. For an interior edge r = ( 7 , 7 ') G Aii, let Afr be the sub¬ 
mesh mad e of t he faces cri,cr 2 of Ai containing r. Let Sk^r) = kerT^ y n 0,-(Zfc) (see 
Dehnition 4.11). The elements of Skin) correspond to splines of Sl{Afr), which are in 
the kernel of Ty^y and with a support in Thus, Skir) C kerT^' for 

any edge r' G Aii, distinct from r. We deduce that any element of Tfc(r) satishes the 
gluing condition along all edges of Ali, and thus corresponds to a spline function in 
Sl{Ai). In other words, we have Tfc(r) C S^(Ai) fl kerTo = /Cfc. The elements of Tfc(r) 
have a support in 7^°'^(r) © 7^°'2(r) and their Taylor coefficients at the end points of r 
vanish. 

Thus the support of the elements of Tfc(r) and Tfc(r') for two distinct edges r, r' do 
not intersect, and their sum is direct. Let Ek = ©reA(i<^fc('r)- 

Let Tk = kerTi fl Sl{Ai) be the set of spline functions, which Taylor expansions 
along all edges vanish. 


Proposition 6.2. 


F'k Ak © Ek 
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Proof. Let / G /C^ and take an interior edge r G Mi. Let cTi,cr 2 be the two faces of 

K. 

Then {fai,fa 2 ) G 5fc(A/'r) Piker Tq =/Cfe(r). By Lemma |4.12| (/^i,/az) = Sr + (^i, A: 2 ) 
with Sr G Tfe(r) and {ki, ^ 2 ) G {Ki, K 2 ). As Sr lifts to a spline G Sl{M), f — Sr is an 
element of Sl{M), which image by the Taylor expansion Tr along the edge r vanishes. 

If r is a bonndary edge of M, which belongs to the face ai, we have a similar 
decomp ositio n fi = Sr + fci with Sr G Tfc(r) and ki G Ki, nsing the convention of 
Remark 4.14 Similarly Sr lifts to a spline G Sl{M), f — Sr is an element of Sl{M) in 


the kernel of Tr. 

Repeating this process for all edges r G Mi, we can constrnct an element / = 
/ — Ylr&Mi ^ A1i,Tt-(/) = 0, i.e. / belongs to kerTi = This 
shows that fCk <Z J^k + constrnction, we have J^k C M and £k = 

®r&Mi£k{,'T) C /Cfc- Considering the snpport of the fnnctions in Jhc and £k, we dednce 
that their snm is direct and eqnal to □ 


6.3. The dimension formula. We can now determine the dimension of S^(M). 
Theorem 6.3. Let s* = max{5(r) | r G Adi}. Then, for k ^ s*, 

dimSliM) = (fc-3)2Fn + |(fc-5)(fc-4)FA 

+ 'Thr&Mi ^T-ik) + 4Tb + 3 Ta — 9Ti + 3To + 


where 


• dr{k) is the dimension of the syzygies of the gluing data along r in degree ^ k, 

• Fa is the number of rectangular faces, F^ is the number of triangular faces, 

• Fi is the number of edges, 

• Fq (resp. Fj^) is the number of (resp. crossing) vertices. 

Proof. By dehnition, we have 


dim 5^ (Ad) = dim.PL + dim/Cfc. 


By Proposition |6.2[ we have 


dim/Cfc = dim Afc + dim Tfc = dim Afc + dimTfc(r) 

r&Mi 

= {k — 3)^Tb + —(fc — 5)(fc — 4 )Ta + ^ } {dr{k) — 9 + Ct-( 7) + Cr'(70) 

tgMi 
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From Proposition 6.1[ we deduce that 


dim5^(A^) 


dim /Cfc + dim "H 

(fc-3)2Fn + i(A;-5)(fc-4)FA 

+ ^ (dr(fc) — 9 + Ct-(7) + Ct-/(7')) 

T=('y,'y')€A 4 i 


+ ^ (F(7) + 3) - 5^ 5^ c,(7) + c+(7) 

7SAI0 ■y&Mo '1&M0 

= Fu{k - 3)2 + Fa^(A; - 5)(A; - 4) + ^ d^{k) - 9Fi 

rSA4i 

+4Fg + 3 Fa 4“ 3-^0 4“ -f+ 

since Er=(7,7')eA4i('^r(7) + Cr'(y)) = E7e^o^r97‘^-^9') and E 7 e^o^(^) = 
3Fa. 


4Fb + 
□ 


As a direct corollary, we obtain the following result: 

Corollary 6.4. If M is a topological surface with gluing data satisfying the compati¬ 
bility Conditions 2.6 2.7 and if all crossing vertices of Ai have 4 edges, then Sl{Ai) 
is an ample space of differentiable functions on Ai for k ^ s*. 


6.4. Basis. We are going now to describe an explicit construction of spline functions 
which form a basis of Sl{Ai). An algorithmic description of the computation of the 
Bernstein coefficients of these basis functions is provided in Appendix A. 

We assume that k is bigger than the separability s* of all edges. 


6.4.1. Basis functions associated to a vertex. Let 7 G A4o be a vertex and cxi,..., ctf 
be the faces of AI 2 adjacent to 7 . We also assume that cxj and cTi_i share the edge 
Ti E Ail and that Ti is not a crossing edge at 7 if such an edge exists. 

To compute the basis functions attached to 7 , we compute first the Taylor coefficients 
of /o-. = p 4- QiUi + qi+iVi + SiUiVi + • • • solutions of the system (26)-(27) and then lift 
these Taylor coefficients to define a spline function with support in (P.y. This leads to 
the following type of basis functions: 

• 1 basis function attached to the value at 7 : p = 1 , = 0 , s* = 0 

• 2 basis functions attached to the derivatives at 7 : p = 0 , [gi, g 2 ] G {[ 1 , 0 ], [ 0 , 1 ]} 
and Sj = 0 if r* is not a crossing edge at 7 and determined by the relations 
(26)-(27) if Ti is a crossing edge at 7 . 

• .^( 7 ) — Xlili ^ri( 7 ) 4 -c+( 7 ) basis functions attached to the free cross derivatives, 
with p = 0, gi = 0 and Sj G (0,1} if r, is not a crossing edge and determined 
by the relations (26)-(27) if r, is a crossing edge at 7 . 


6.4.2. Basis functions associated to an edge. Let r be an edge of Aii shared by two 
faces (Ji, (72 with vertices 7 , 7 '. Let us assume that the coordinates of these points in 
the face ai are 7 = ( 0 , 0 ) and 7 ' = ( 1 , 0 ). 

The elements of are the image by 0,- of the elements of of the form 

P [Ai,Bi,Ci]+Q [A 2 , 82 , 02 ] with degree deg(P) ^ k-fi^-mr, deg(Q) ^ k-u^-mr 
which are in the kernel of and T^i. 
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From the relation (23), we deduce that P(0) = 0, Q(0) = 0. That is, P and Q are 
divisible by u. 

• If Cr(7) = 0 , i.e. 7 is not a crossing vertex, we have i?i( 0 )C' 2 ( 0 ) —-B 2 ( 0 )C'i( 0 ) = 

a( 0 ) 7 ^ 0 and the relation (|^ implies that P'{0) = 0 , Q'( 0 ) = 0 . That is 

P = u^P, Q = u^Q. 

• If Ct-( 7 ) = 1, then the kernel of is generated by polynomials such that -P(O) = 

0, Q(0) = 0, P'(0) = AC 2 ( 0 ), Q'(0) = -AC'i(O). That is P = u{\C2{0)+uP), 

Q = u (-ACi(O) +uQ). 

That is 

P = u {x c^( 7 ) ( 72 ( 0 ) +uP^ ,Q = u ^-A Cr( 7 ) Ci(0) +uC^ . 

By symmetry at 7 ', we see that P and Q are of the form: 


P = M (1 - m) ^A c^( 7 )C' 2 ( 0 ) (1 - m) + A'c^(7') 6*2(1) M + M (1 - m)-P j , 

Q = —u (1 - m) ^A c^(7) 61(0) (1 - m) + A'c^(7') 61(1) M + M (1 - m) , 


with A, A' G M, deg(P) — /i — m — 4, deg((5) ^k — u — m — A. 

We construct a basis of Tfc(r) by taking the image by 0 ,- of a maximal set of linearly 
independent elements of this form (see Section 4.3). This yields (iT-(/c) — 9 + Ct-( 7 ) + Ct-( 7 ') 
spline basis functions. 


6 . 4 . 3 . Basis functions associated to a face. Finally, we dehne the basis functions at¬ 
tached to a face a G M.2 as the 2 -interior Bernstein basis functions in degree ^ k. 
There are [k — 3 )^ such basis spline functions for a rectangular face and (^2^) for a 
triangular face. 


7. Examples 

7 . 1 . Splines on flat triangular tilings. We consider a subdivision of a planar do¬ 
main n C into a partition of triangles and the topological surface M. induced by 
this subdivision. 

For two faces cri,cr2 G AI2, which share an edge r G M.\ at a vertex 7, there is 
a linear map 00-2,0-1, which transforms the variables (Mi,ni) attached to cxi into the 
variables {u2,V2) attached to (T2. 

With 7 = ( 0 , 0 ) and vi = M2, the transition map 0o-2,o-i is given by 


U2 


■ 0 

b ■ 


Ml 

. ^2 . 


1 

a 


. . 


where a, 6 G M and 6 7^ 0 . We choose these constant transition maps to define the 
space of splines iS^(Al). The gluing conditions along the edges correspond then to 
conditions for the polynomials expressed in the same coordinate system. In this case, 
5 ^(A 1 ) is the vector space of piecewise polynomial functions on A 4 , which are on 
n, that is, the classical C'^-spline functions on hi. 

If a = 0 , the edge of (J2 at 7 distinct from r is aligned with the edge of ui at 7 
distinct from r. The vertex 7 is a crossing vertex (c+(7) = 1 ; all the coefficients a 
in the transition maps around 7 vanish) if there are 4 edges at 7, which are pair-wise 
aligned. 











26 


B. MOURRAIN, R. VIDUNAS, AND N. VILLAMIZAR 


As for any interior edge r G Mi the transition map is constant, we have n,- = 0, 

jUt- = 0, = 0, = 1, ^ 5 and = 2k. For the bonndary edges, we have 

dr{k) = 2fc + 2 (see Remark 4.14). 

We dednce from Theorem 6.3| that for fc ^ 5, we have 

dim = 2 — 4)FA + 2k F° + (2k + 2) + 3FA — 9 Fi + 3 Fq + F+ 

= — (fc + 2) (F + 1) FA — G (k — 2) FA + (2k — 9) F° + (2k — 7) F’{’ + 3 Fg + Fl|_ 

where F° (resp. Fj^) is the nnmber of interior (resp. bonndary) edges. Using the 
relations 3 Fa = 2F° + Fj’ (connting the edges per triangle, we connt twice the interior 
edges shared by two triangles and once the bonndary edges), Ff = Fq and Fq = Fq+Fq 
where Fq (resp. Fq) is the nnmber of interior (resp. bonndary) vertices, we obtain 

dimSl = ^(F + 2)(F + 1)Fa-( 2F + l)F° + 3Fo" + F+. 

This coincides with the dimension formnla of piecewise polynomials of degree F ^ 5 
on a triangnlar planar mesh, given in HZI. Here F+ connts the nnmber of crossing 
vertices, also called singnlar vertices in HU. 

The basis fnnctions constrncted as in Section [R4l are as follows: 

• For each vertex 7, there are 3 basis fnnctions associated to the evaluation and 

derivatives in x,y at 7. There are ^(7) — u(7) + ^+(7) basis functions 

associated to the free cross-derivatives on the triangles containing 7. 

• For each interior edge r = (7,7'), there are 2F—9-|-CT-(7)-l-Cr(7') basis functions 
associated toF-|-l — 6 = F — 5 free interior Bernstein coefficients ..., &fc- 3 ,o 
on the edge, F — 4 free interior Bernstein coefficients 62,1, • ■ ■, Ffc-2,1 on one 
triangle a which contains r and 62,0 (resp. Ffc-2,0) if t is a crossing vertex at 7 
(resp. 7'). 

• For each boundary edge r', there are 2k —7 basis functions associated to F — 3 
free interior Bernstein coefficients 62,05 ■ ■ ■ 5 bk- 2,0 on the edge r', and F — 4 free 
interior Bernstein coefficients 62,15 • • • 5 6^-2,1 on the triangle a which contains 
t' . 

• For each triangle a, there are (^2^) basis functions associated to the interior 
Bernstein coefficients bij with 2^i,j^k — 2 and 0 ^ i + j ^ k. 

This basis description involves the Bernstein coefficients of polynomial on the triangles. 
The basis differs from the nodal basis proposed in HI]. From the listed Bernstein 
coefficients, we can however recover the nodal basis of im. dual to the evaluation and 
derivatives at the vertices and at interior points of the edges and the triangles. 

7.2. A round corner. We consider a mesh M composed of 3 rectangles cti,(T2,ct3 
glued around an interior vertex 7, along the 3 interior edges ti,T2,T3. There are 6 
boundary edges and 6 boundary vertices. 

We take symmetric gluing data at 7 and at the crossing boundary vertices Si. The 
transition map across the interior edge Xj is given by the polynomials: [a, 6, c] = [(u — 
1), — 1,1] where 7 is the end point with u = 0 and Si is the end point with u = 1. The 
generating syzygies are 


5l = [0,l,l],^2=[l5W5l]. 
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£3 



Figure 5 . Smooth corner. 


For the interior edges r*, we have n = l, m = 0, /i = 0, i/ = l and = k + l + k = 

2k + 1 . For the boundary edges r', we have n = 0, m = 0, u = 0, i/ = 0 and 
drik) = 2{k + l). 

As a(0) = —1 (resp. a(l) = 0), 7 is not a crossing vertex {Cni'j) = 0) and 6 i is a 
crossing vertex of r* (cT--( 5 i) = 1). 

We check that the separability of all the interior edges is 4. For fc = 4, the dimension 
of SUM) is 

3x(4-3)^ + 3x(2x4 + l) + 6x(2x4 + 2)+4x3-9x9 + 3x7 + 6 = 48. 


The basis functions are constructed as in Section 6 . 4 , using the algorithms of Appendix 
A. 


• The number of basis functions attached to 7 is 6 = 1 + 2 + 3 . 

— The basis function associated to the value at 7 is 

[^0,0 + ^1,0 + ^0,1 + ^1,1; bofi + &i,o + ^0,1 + ^1,1) ^0,0 + ^1,0 + ^0,1 + ^1,1] 


— The two basis functions associated to the derivatives at 7 are 

[I ^1,0 + i ^1,1 + ^ ^2,0 + 12 ^2,1 + I ^1,25 

— \ bo,I — \ bi,i — I 62,1 — ^ bo,2 — bi, 2 , 

I bo,i - i bi,o - ^ 62,0 - ^ b2,i + ^ bo,2 + I 61,2] 

[i ^0,1 + i ^1,1 + I ^2,1 + ^ bo,2 + &1,2, 

— \ bo,i + I bi,o + ^ &2,o + ^ ^2,1 — bo ,2 — I bi, 2 , 

bi,o - I bi,i - b 2 ,o - ^ &2,i - I ^1,2] 

— The three basis functions associated to the cross derivatives at 7 are 

[ik ^1.1 “ i^ ^2,0 - ^ &2,i - ^ bo,2 - I bi,2, 
b2,o - 17 &2,i, -]7 bo,2 - I bi,2] 

[“17 bo,2 — I bi,2, 

ik ^1-1 ~ 17 ^2,0 - ^ b2,i - 17 bo,2 - | &17, “17 &2,o - 17 ^2,1] 

[“17 ^2,0 - 17 ^2,1, 

“17 bo,2 - I bi,2, 17 61,1 - 17 62,0 - ^ b2,i - 17 &o,2 - | ^1,2] 
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• The number of basis functions attached tohj is 4 = 1 + 2 + 2 — 1. Here are 
the 4 basis functions associated to 

[^3,0 + ^3,1 + ^4,0 + ^4,1,0, 6 o ,3 + ^1,3 + ^0,4 + ^1,4], 

[^3,0 + ^3,1,0, 60,3 + ^1,3], 

[^3,1 + ^4,1; 0, —&1,3 — 

[^3,1, 0, —fel,3]- 

The basis functions associated to the other boundary points 62 , S 3 are obtained 
by cyclic permutation. 

• The number of basis functions attached to the remaining boundary points is 
4 = 1 + 2 + 1. For ci, the 4 basis functions are 


[^3,3 + ^3,4 + ^4,3 + ^4,4, 0, 0], [&3,3 + ^3, 0, 0], [&3,3 + ^3,4) 0) 0]) [^3,3, 0, 0] 


The basis functions associated to the other boundary points are obtained by 
cyclic permutation. 

• The number of basis functions attached to the edge r* is 2x4 — 7 = 1. For 
the edge ri, it is 

[^ 2 , 1 , 0 , —& 1 , 2 ]- 

The basis functions associated to the other interior edges are obtained by cyclic 
permutation. 

• The number of basis functions attached to the boundary edges is 2(4 — 3) = 2. 
For the boundary edge (ei, 5i) of cxi, the two basis functions are 

[^3,2, 0, 0], [64^2, 0, 0]. 

• The number of basis functions attached to a face dj is (4 — 3)^ = 1. The basis 
function associated to ai is 


[62,2, 0, 0] 


and the two other ones are obtained by cyclic permutation. 


7.3. A pruned octahedron. We consider a mesh Ai with 6 triangular faces AEF, 
CEF, ABE, BCE, ADF, CDF and one rectangular face ABCD, depicted in Figure]^ as 
the Schlegel diagram of a convex polyhedron in It is an octahedron where an edge 
BD is removed and two triangular faces are merged into a rectangular face (see [22] for 
the complete octahedron, which involves only triangular faces). 



Figure 6. A pruned octahedron 
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We are going to use the following notation for the variables on the faces of M. For 


X,Y G two vertices dehning an edge XY of a face cr, let u 


-)■ 


be the 


linear function with u\{X) = 0, u\{Y) = 1 and u\{Z) = 0 for all the points Z on 
the edge of a through X and distinct from XY. We will use these linear functions 
tif, etc., as variables on the different faces. As the restriction on a share edge 
XY of the two functions dehned on the faces adjacent to XY coincide, there is no 
ambiguity in evaluating these linear functions on Ai. For a triangular face XYZ, we 
have + Uy + = 1. For a rectangular face XYZW, we have and 

= Uy. We denote by dxy the derivative with respect to the variable u^. It is such 
that dxY{u^) = 1- On a triangular face XYZ, we have dxy + dyz + dzx = 0. 

We use a symmetric gluing at all the vertices and therefore the vertices A, C, E, F 
are crossing vertices. Let us describe how we construct the gluing data on the edges 
by interpolation at the vertices, in a smaller degree than the degree associated to the 


gluing (14) proposed in na. 

In terms of differentials (see relation (|^), the symmetric gluing at the vertices trans¬ 
lates as 



Bea + Bec = 

0 , 

Beb 

+ Bee — 

0 

at 

E, 


Bea + Bfc = 

0 , 

Bed 

+ Bee = 

0 

at 

F, 

(28) 

Bab + Baf = 

0 , 

Bae 

+ Bad = 

0 

at 

A, 


Bcb + Bcf = 

0 , 

Bce 

+ Bcd = 

0 

at 

c. 

and around the vertices of order S 

) . 







Bba + Bbc + Bbe — 0 

at B, 


Bda + Bbc 

+ Bdf 

= 


For gluing the triangles EFA and EEC along EE, we interpolate the following relations 
between the derivatives: 

Sea + dEc = 0 at E, Oea + Qec = 2 Bee at F, 

where the second expression is Of a + Bfc = 0 rewritten using Of a = Bea — Bee, Bfc = 
Bec — Bee- We choose the linear interpolation 

Be A + Bec = 2^^ Bee- 

Thereby we have = 2 m;|, ^ef = —1 and the gluing data for the edge EE is 

[ 2 n^,-l,l]. 

For the edge EB between the triangles EBA and EBC, we interpolate the following 
relations: 

Bea + Bec = ^ at E, Bea + Bec = 3 Beb at B, 

where the latter relation is Bba+Bbc+Bbe = 0 rewritten using Bba = Bea—Beb, Bbc = 
Bec — Beb, Bbe = —Beb- Additionally, we have to take into account the compatibility 
conditions (|^-(10) since is a crossing vertex. It translates as Beb{<^eb) = Bef{<^ef) 
and BEci'^Ec) = Bea{.^ea) at the vertex E. This leads to the following gluing data on 
the edge EB: 

EB: [2m| + (m|)2,-1,1]. 

Similarly, the gluing data of the edge FD is 

FD:[2u^ + {u^f,-l,l]. 
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The edges EA, EC, EA, EC connect cross vertices just as EE and yield linear gluing 
data 


EF: [ 24 ,-1,1], EC-. [2u%,-l,l], 
FA: [2n;^,-l,l], FC: [2n^,-1,1] 


We check that the compatibility conditions ([^-(10) are satished across EA,EC and 
FA, FC. The gluing data along AB, AD, CB, CD looks the same: 


AB: 

CB: 




AD: 

CD: 


[2ng,-l,l]. 


We have linear gluing data everywhere except on the edges EB and FD. Let us analyze 
the syzygies associated this data. 

• For the edges EB and FD with one crossing vertex, the gluing data is of the 

form [2u + —1,1]. We have n = 2 and m = 1 since the edge is connecting 

two triangles, /x = 0 and v = 2 and d{k) = 2k — 2. The /i-basis is [0,1,1], 
[—1, —2u — V?, 0]. The separability is achieved in degree /c ^ 6 and not 5 as it 
could be expected (d(5) ^ 8 ). 

• For the edges EA, EA, FA, FC, FF between triangular faces, with two crossing 
vertices, the linear gluing data is of the form [2u, 1, —1]. We have n = 1 and 
m = l, fi = 0 , 1^ = 1 and d{k) = 2k — 1. The /i-basis is [0,1,1], [—1, —2u, 0]. 
The separability is achieved in degree k ^ 4. 

• For the edges AB, AD, CB, CD between a triangular face and a rectangular 
face, with one crossing vertices, the linear gluing data is of the form [ 2 m, 1 , — 1 ], 
We have n = 1 and m = 0, /i = 1 since the degree of the homogeneization 
[da, dfo, dj (see Dehnition 4.1) is [3,2,1] or [3,1,2], v = 1 and d{k) = 2k. The 
/i-basis is [0,1,1], [—1,—2 m, 0]. The separability is also achieved in degree 
k>4. 


Now we count how many splines do we have in degree k ^ Q: 

• For the four crossing vertices A, C,E,F we have 1 -|- 2 -|- 1 = 4 dimensions and 
1 -|- 2 -|- 3 = 6 dimensions for B and for D. In total we have 4 ■ 4 -|- 2 ■ 6 = 28 
degrees of freedom around the vertices of Ai. 

• For the edges EB and FD, we have 2(fc — 2) — 8 = 2fc — 12 dimensions. For 
the edges EA, EA, FA, FC, EF, we have 2k — 1 — 7 = 2k — S dimensions. For 
the edges AB, AD, CB, CD, we have 2k — 8. 

• For the 6 triangular faces, we have (^ 2 ^) dimensions and for the rectangular 
face {k — 3)^. 

The dimension formula in degree fc ^ 6 is then 


-4 +11-2(^-4)+ 6 


k-4 
2 


+ {k-3f 


{2k-3f + k-4. 


For k = Q, the dimension is 83. It turns out that this formula also holds for degree 
k = 4, k = 5. 

The construction of basis functions can be done as described in Section 16.41 Let us 
give the basis functions associated to the value and hrst derivatives at the point A. 
Here are the Bernstein coefficients in degree 4 of the basis function for the value at A 
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with the vertex A represented in the center and the edges represented by horizontal 
and vertical central lines (in bold): 


0 

0 


0 

0 

-1 

0 


0 

1 

1 

1 

-1 

0 


0 

1 

2 

1 

1 

1 

0 

0 


0 

1 

1 

1 

0 




E 


0 

1 0 

0 

B 

A 

F 


D 


This gives the following specializations to the polygons ABE, AEF, AED, ABCD (re¬ 
spectively), selectively de-homogenized: 

(l + 3«I)(l + 2nf-n^), 

{'^eY (1 + + 2 m® -|- Qu^u^), 

{upY (1 + 3m5)(1 + 2m® — M®), 

{u^u^f + ?>u^) - 2Au^u^{u%u^ + u^u^)), 

and 0 on the other faces. The basis function associated to the hrst derivative in one of 
the directions at the cross vertices is: 


(29) 


0 

0 


0 

0 

7_ 

24 

0 


0 

1 

3 

0 

16 
_ 1 
6 
0 


0 

1 

6 

1 

4 

0 
_ 1 
4 

0 

0 


0 
1 
3 
0 
_ 1 
3 

0 


0 

0 

0 


0 


The non-zero specializations to ABE, AEE, AED, ABCD are, respectively: 

( mb )^ m ®(1 + 3 m ®), ( m ^)^ m ®(1 + 3 m ®), -{upf {up + Au^), 

{uiu^ fu^ + 3u^) + 7u^uiu^) . 

The basis function for the derivative in the other direction is obtained by a mirror 
image of (29). 

The basis function corresponding to the cross derivatives is realized by 

0 


0 

0 



0 

0 

0 


0 

1 

12 

0 

1 

12 

0 

0 

0 

0 

0 

0 

1 

1 

0 

1 

0 

24 

16 

12 

0 

1 

24 

0 

0 



0 

0 




0 


In this spline, we could modify the 0 entry next to two A entries to A^ so to lower the 
degree of the specialization to the rectangle. After the modihcation, the 4 non-zero 
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specializations would be 




The local splines around C look the same. The local splines around other vertices 
involve the edges EB and FD, and we would need degree 6 splines. 
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Appendix A. Algorithms for the basis construction 

Our input data is the topological surface Ai and the gluing data. For each edge r 
oi M., we are given the /i^-basis of Z^. 

SI = [A[, Bl, or], ^2^ = [A^, Bl, OJ]. 

The rational map is then described by ar = —,br = ^ with Ur = BIC^ — B^Cf, 
K = AlCl - AlCl, and - AlBf. 

The spline basis functions / = {fa) are represented on each face a by their coefficients 
in the Bernstein basis of the face in degree k: 


Let ea{k) = if a is a rectangular face and ea{k) = fc (fc — 1) if cr is a triangular face. 


A.l. Vertex basis functions. Let 7 be a vertex of M. shared by the faces ai,... ,(jp 
and such that at and cxi+i share the edge Tj+i. We compute the Bernstein coefficients 
= [co,oAi oAo 1 Ai,i 5 • • •] of fhe basis functions attached to a vertex at, using the 
equations (26), (27) and the relation between the Bernstein coefficients and the Taylor 
coefficients of the function at (0,0), see (16), (17). 

If Cg Q corresponds to the point 7, with coordinates (0, 0) in the face cxj and /o-^ = 
p + QiUi + qi+iVi + SiUiVi + ■ ■ • are c*, we use the relations p = Cg g, g* = kc\Q, = 

Si = efc(ci_i- 4 g- 4 i). 
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Basis function for the value at vertex 7 


for i in [1,F] do 

let := 1 , 4 0 := 1 , 4 ,i := 1 , := 1 and 4 ^ := 0 for 

(/,m) ^{(0,0),(1,0),(0,1),(1,1)}; 

end 

Basis functions for the derivatives at vertex 7 

for [4o>Co,i] in {[1,0], [0,1]} do 
for i in [2,F] do 



’ 4,0 ’ 

_ I 0 1 u 


. 4,1 _ 

1, 1J„(0) a,(0) ) [ 


v-i 

^ 1,0 

V-l 

^ 0,1 


end 

if all edges Ti are erossing edges at 7 then 
let c\ i := 0 ; 

end 

for i in [2,F] do 

if Ti is a erossing edge at 7 then 

^1,1 ~ ^1,0 “T Cq i 

(47i^ - 47o^ - cl~i) 

+k 440) 477 + k b7(0) c \ J ) 

else 

I 4,1 := 0; 

end 
end 

for i in [1,F’] do 

lift(Ci_i,Ci, Ti) 

end 
end 

Basis functions for the cross derivatives around 7 
for i in [1,F] do 

let := 0,4 o := 0 ,c 4 := 0 ; 

end 

if all edges Ti are erossing edges at 7 then 
let cj 1 := 1 ; 
for i in [2,F] do 

4,1 = £^^.(0) 474 

end 

for i in [1,L] do 

lift(ci,Ci_i, Ti); 

end 


else 


for i in [1,F] sueh that Tj is not a erossing edge at 7 do 
let 4 1 = 1 and c 4 = 0 for / 4 J j 
for i in [1,F’] do 

if Ti is a erossing edge then 


4,1 = S)^’-4o)4,i' 
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The function lift(cj, Cj_i, Ti) used in the algorithm consists in computing the coefficient 
of a spline function with support along the edge Tj, from its hrst Taylor coefficients on 
the faces <7*. 

LIFT(Ci,Ci_i, Ti) 
for i in [ 1 ,F] do 
solve the systems: 

r kci, 1 ^ r Ai(o) ^2(0) 1 r ■ 

L \ [ 5 i( 0 ) 52(0) \[ Qo _ 

and 

{ci i — Cifl ~ Cp 1 ) _ -Bi(O) -^2(0) Po 

(ci 1 — Q — Cq _ C'i(O) ^*2(0) . . 9 'o _ 

_ r Bi{ 0 ) ^2(0) 1 r p\ 

[ C'i(O) ^ 2 ( 0 ) \[q{ _■ 

compute P* := »n(l — Suf + 2 uf) + p](u — 2u‘f + uf), 

Q" ■= go(l - + 2 Mf) + q\{u - 2 uj + u^); 

compute the image {pi, Pi) of P^S\ + Q'^S\ by ©t; and update the coefficients of 

^i—li ^ii 

end 


As Ai( 0 )P 2 ( 0 ) — A 2 ( 0 )Pi( 0 ) = c(0) 7 ^ 0, the hrst system has a unique solution. When 
Pi( 0 )C 2 ( 0 ) — P 2 ( 0 )Ci( 0 ) = a( 0 ) 7 ^ 0 (i.e. when r* is not a crossing edge at 7 ), the 
second system has a unique solution. When a(0) = 0 (i.e. when r* is a crossing edge 
at 7 ), the second system is degenerate, but it still has a (least square) solution. 

The polynomials P* (resp. Q*) are constructed so that P*(0) = Po,P* (0) = p\, 
P*(l) = 0,P''(1) = 0 (resp. g'(0) = qlQ^'{0) = q{, Q^{1) = = 0). 


By construction, the Taylor expansions of their image by 0^-^ vanish at 7 ' and coincide 
with [co 7 o\ c; 7 o\ 

at 7 respectively on (Tj_i and Uj. 


*^1,0 *^0,1)]) [*^0,0’^i,o> *^0,0 ^1,0 '^0,1). 
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A.2. Edge basis functions. 

Basis functions for the edge r 

Inpnt: the /i-basis of the syzygy module Z(t); 

if Ct-( 7) = 1 then 

compute the image by 0^ of u(l — (<^ 2 ( 0 ) ^ 1 ] 

end 

if Cri'j') = 1 then 

compute the image by Qr of m ^(1 — u) (CJ(1) [A[, Ci] 

end 

Let A = u^{l — u)^; 
for i in [0,k — /x — m — 4/ do 

compute the image by 0,- of u*A [A[, i?[, Cf]. 

end 

for i in [0,k — u — m — A] do 

compute the image by 0,- of u*A [AJ, -BJ, CJ]. 

end 

A.3. Face basis functions. 

Basis functions for the face a 

for 2 ^ i ^ k — 2, 2 ^ j ^ k — 2 (and i+j^k — 2ifaisa triangle) do 
let Cij := 1 and Ci>ji = 0 for i or j' 7 ^ j. 

end 
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